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Abstract. We numerically investigate the transport properties of interacting spinless elec- 
trons in disordered systems. We use an efhcient method which is based on the diagonalization 
of the Hamiltonian in the subspace of the many-particle Hilbert space which is spanned by 
the low-energy Slater states. Low-energy properties can be calculated with an accuracy com- 
parable to that of exact diagonalization but for larger system sizes. The method works well 
in the entire parameter space, and it can handle long-range as well as short-range interac- 
tions. Using this method we calculate the combined effect of disorder and interactions on the 
Kubo-Greenwood conductance and on the sensitivity of the ground state energy to a twist 
in the boundary conditions. We find that the influence of the interactions on the transport 
properties is opposite for large and small disorder. In the strongly localized regime (small 
kinetic energy, large disorder) interactions increase the transport whereas for weak disorder 
(large kinetic energy) interactions decrease the transport. 
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1 Introduction 

The influence of electron-electron interactions on Anderson localization has reattracted 
a lot of attention, in particular because recent experimental results show a metal- 
insulator transition (MIT) in the two-dimensional electron gas in Si-MOSFETs [Q. 
This is in conflict with the theory of Anderson localization for non-interacting elec- 
trons. Since the electron density in the Si-MOSFETs is very low which makes the 
electron-electron interaction particularly important, it is generally assumed that some 
type of interaction effect is responsible for the MIT. A complete understanding has, 
however, not yet been obtained. There have been attempts to explain the experiments 
based on the perturbative renormalization group H, non-perturbative effects [pi, or 
the transition being a superconductor-insulator transition rather than a MIT HT 

We have investigated the influence of disorder and interactions on the transport 
of spinless electrons numerically. Within an effective single-particle theory at the 
Hartree-Fock level interactions always have a localizing effect g . In order to go beyond 
the Hartree-Fock approximation we have developed Iq, m an efficient method, the 
Hartree-Fock based diagonalization (HFD) which is related to the quantum-chemical 
configuration interaction approach. We have used this method to study the inffuence 
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of interactions on the conductance in one ||], two |^, and three Q dimensions. We 
found a delocahzing tendency of the interactions for strong disorder but a locaUzing 
one for weak disorder. Similar resuhs have been obtained by means of the density- 
matrix renormahzation group llG] in one dimension and exact diagonahzation in two 
dimensions [Ol. 

In this paper we compare the results for the Kubo-Greenwood conductance to 
those for the phase sensitivity, i.e., the reaction of the ground state energy to a twist 
in the boundary conditions. We only summarize the basic findings, the details will be 
published elsewhere. 



2 Quantum Coulomb glass and Hartree-Fock based diagonalization 

The generic model for spinless interacting disordered electrons is the quantum Coulomb 
glass P, 113]. It is defined on a regular hypercubic lattice with g — L''' {d is the 
spatial dimensionality) sites occupied by TV = Kg electrons {0 < K <1). To ensure 
charge neutrality each lattice site carries a compensating positive charge of Ke. The 
Hamiltonian is given by 

H = -tY,i4cj + c]ci) + ^ if.n, + - ^(n, - X)(n, - K)U,, (1) 

(y) « i#J 

where cj and Ci are the electron creation and annihilation operators at site i, respec- 
tively, and {if) denotes all pairs of nearest-neighbor sites, t is the strength of the 
hopping term, i.e., the kinetic energy, and rn is the occupation number of site i. We 
parametrize the Coulomb interaction C/^ — e^ jr^ by its value U between nearest- 
neighbor sites. The random potential values tpi are chosen independently from a box 
distribution of width 2Wo and zero mean. The boundary conditions are periodic 
with an additional Bloch phase in one of the dimensions. The Coulomb interaction is 
treated in the minimum image convention. 

A numerically exact solution of the quantum Coulomb glass requires the diagonal- 
ization of a matrix whose dimension increases exponentially with system size. This 
severely limits the possible sample sizes. In order to overcome this problem we have 
developed the HFD method. The basic idea is to work in a truncated Hilbert space 
consisting of the corresponding Hartree-Fock (Slater) ground state and the low-lying 
excited Slater states. For each disorder configuration three steps have to be performed: 
(i) find the Hartree-Fock solution of the problem, (ii) determine the B Slater states 
with the lowest energies, and (iii) calculate and diagonalize the Hamilton matrix in 
the subspace spanned by these states. The number B of new basis states determines 
the quality of the approximation, reasonable values have to be found empirically. 



3 Kubo-Greenwood conductance and phase sensitivity 

In order to characterize the transport behavior we have studied two quantities, the 
Kubo-Greenwood conductance and the phase sensitivity of the ground-state energy. 
The conductance |1^, |lj| can be obtained from linear-response theory. It is essentially 
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Fig. 1 Logarithmically averaged (400 samples) d.c. conductance G(0) for systems of 4 x 4 
lattice sites and 8 electrons for different U and t. The disorder strength is fixed to Wo = 1, 
the broadening is 7 = 0.05, and the HFD basis size is S = 300. 



given by the current-current correlation function of the ground state. For an isolated 
finite system the conductance as a function of frequency consists of a finite number of 
5 peaks at the excitation energies of the system. In order to extrapolate to the d.c. 
conductance at zero frequency an inhomogeneous broadening 7 is assumed in our cal- 
culations which mimics the coupling to leads and contacts. This broadening introduces 
an additional empirical parameter. We have calculated the conductance for systems 
with up to 100 lattice sites. In Fig. |l| we show a typical result (for a comparatively 
small system which also allows the calculation of the phase sensitivity). For weak 
disorder (large kinetic energy t) the interactions always reduce the conductance while 
for strong disorder (small t) moderate interactions significantly increase the conduc- 
tance. Sufficiently strong interactions always strongly suppress the conductance. This 
is the precursor of a Wigner crystal or Wigner glass. The behavior of the conductance 
can be attributed to the competition of two effects: First, the interactions destroy 
the phase of the electrons and thus the interference necessary for localization. This 
is particularly effective if the localization length is small to begin with. Second, the 
interactions introduce an additional source of randomness which tends to increase the 
localization. 

It is well known ||l^ that in disordered systems there is a close relation between 
the charge stiffness (i.e., the 2nd derivative of the ground state energy with respect 
to a magnetic fiux) and the conductance. A numerically simpler but related quantity 
is the phase sensitivity of the ground state. It is given by the difference between the 
ground-state energies for periodic and antiperiodic boundary conditions [|l6| . In Fig. Q 
we show numerical results for the phase sensitivity of the same systems as considered 
in Fig. |l|. For strong disorder (small t) the phase sensitivity shows the same qualitative 
behavior as the conductance. However, the interaction-induced enhancement of the 
phase sensitivity for moderate interactions is weaker than the enhancement of the 
conductance. For small disorder (large t) a crossover to a localizing effect of the 
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Fig. 2 Logarithmically averaged (400 samples) phase sensitivity Lx AE. The parameters 
are as in Fig. |l| The basis size within the HFD method was B = 500. 

interactions starts, but to really reach the localizing regime seems to require still larger 
t. The differences between the results for the Kubo-Greenwood conductance and the 
phase sensitivity are not completely understood so far. However, a likely reason are 
the ambiguities in describing the contacts and in taking the zero- frequency limit. 

In summary, we have studied the influence of electron-electron interactions on 
Anderson localization for spinless electrons in two dimensions. For strong disorder 
moderate interactions significantly enhance the transport. 

This work was supported in part by the DFG (SFB 393/C2). T.V. thanks the Aspen Center 
for Physics for hospitality during the completion of this paper. 



References 



[8' 

[to 

11 
12' 
13 
14 
15; 

[16] 



S. V. Kravchenko, D. Simonian, and M. P. Sarachik, Phys. Rev. Lett. 77 (1996) 4938 

C. Castellani, C. DiCastro, and P. A. Lee, Phys. Rev. B 57 (1998) R9381. 
S. Chakravarty, S. Kivelson, C. Nayak, and K. Volker, |cond-mat/9805383| 

D. R elitz and T. R,. Kirk patrick. Phys. Rev. B 58 (1998) 8214; F. C. Zhang and T. M. 
Rice, |cond-mat/9708050t P. Phillips ei al, Nature 395 (1998) 253 
F. Epperiem, ' " ' "' " ^ -- ' n 
T. Vojta, " ■ _ , 
T. Voita, F. Epperlein. and M. Schreiber, accepted for publication in Comp. Phys. 



vi. scnreioer, and T. Vojta, Phys. Rev. B 56 (1997) 5890 
F. Epperlein, and M. Schreiber, Phys. Rev. Lett. 81 (1998) 4212 
Vojta, F. Epperlein. and M. Schreiber, accepted for publication in Comp. 
Commun., hond-mat/980917l| 
M. Schreib er, F. Epperlein, an d T. Vojta, Physica A 266 (1999) 443 
T. Vojta and F. Epperlein, Ann. Phys. (Leipzig) 7 (1998) 493 

P. Schmitteckert, II. A. Jalabert, D. Weinmann, and J.-L. Pichard, Phys. Rev. Lett. 81 
(1998) 2308 



;ond-mat/9904096 



G. Benenti, X. Waintal, and J.-L. Pichard, 

J. Talamantes, M. Pollak, and L. Elam, Europhys. Lett. 'Sb (l996) 511 

K. Kubo, J. Phys. Soc. Jpn. 12 (1957) 576 

D. A. Greenwood, Proc. Phys. Soc. 71 (1958) 585 

see, e.g., G. Montambaux in Quantum Fluctuations, Les Houches, Session LXIII, edited 

by E. Giacobino, S. Reynaud, and J. Zinn- Justin, Elsevier, Amsterdam 1995 

Care has to be taken if (avoided) level crossings occur as a function of the boundary 

phase shift. 



